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Abstract 



The formalism of nonequilibrium perturbation theory was constructed by 
Schwinger and Keldysh and then was developed with the diagrammatical 
technique by Lifshitz and Pitaevskii. Until now there has been widespread 
application to various researches in physics, condensed matter, plasmas, 
atoms and molecules, nuclear matter etc.. In spite of this, the formalism 
has not been established as perturbation theory. For example, there is no 
perturbative method to derive arbitrary self-energy properly. In addition, 
the connection with other formalism, e.g., the Matsubara imaginary-time 
perturbative formalism is uncertain. Although there must be the relation- 
ship between self-energies in the perturbative formalism, such basic prob- 
lems remain to be solved. The solution is given by the present work. The 
real-time perturbative expansion is performed on the basis of the adiabatic 
theorem. As the results, the requirements of self-energies as functions in time 
are demonstrated and the formulated self-energies meet the known relations. 
Besides, it gives exact agreement with functions derived by perturbative ex- 
pansion in imaginary-time and analytical continuity. As a consequence, it 
implies that the present formalism can be generalized. 

Next, using the formulated self-energies, the behavior of the Kondo reso- 
nance is investigated for nonequilibrium states caused by bias voltage. As nu- 
merical results, the Kondo peak disappears when voltage exceeds the Kondo 
temperatures; it is supported by experiments for two terminal systems. Over 
ten years, it has been being waited in expectation that the Kondo peak splits 
owing to bias voltage as a candidate for two channel Kondo effect. Never- 
theless, it has not been observed in two terminal systems by experiments. 
Here, it is discussed why the Kondo peak splitting may not arise in normal 
two terminal systems. 
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Chapter 1 



Introduction 



1.1 Nonequilibrium Perturbative Formalism 
( Schwinger-Keldysh Formalism ) 



The basic idea on the nonequihbrium perturbation theory was proposed 
by Schwinger in 1961. [1] That included the essentials for the nonequilibrium 
perturbation theory. The idea of time reversal is the basis of the theory. The 
real time-contour has the positive and reverse time directions, so that it starts 
and ends at t = — oo by way of t = oo. In addition, the matrix form in the 
nonequilibrium Green's functions is written after the real time-contour. After 
that, there are two main developments in the nonequilibrium perturbation 
theory. One is the expansion of the formalism by means of the equation of 
motion in the Green's function by Kadanoff and Baym[2]; there their time- 
contour includes imaginary-time path. Another is that the formalism was 
extended as the frame of the nonequilibrium perturbation theory by Keldysh 
in 1965. [3] The formalism is constructed using density matrix and S-matrix 
after the time-contour and the Dyson's equation in the matrix form of the 
nonequilibrium Green's functions on the basis of the idea of Schwinger. The 
Dyson's equation for the Keldysh formalism is given by 



G 
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Then the Keldysh formahsm was studied in further detail and generahzed 
as the formahsm of the nonequihbrium Green's functions and the perturba- 
tive method with the help of the diagrammatic technique by Lifshitz and 
Pitaevskii. [4] 

The main progress in the formalism of the nonequihbrium Green's func- 
tions and the nonequihbrium perturbation theory is summarized with the 
aid of chart as follows: 



Basic idea on the nonequihbrium perturbation theory 
( the real time-contour and the matrix form ) was proposed 
by Schwinger in 1961 



Formalism was constructed 
after the time-contour 
and the matrix form 
by Keldysh in 1965 



Formalism was expanded 
using the equation of motion 
by Kadanoff and Baym 



(1961-62) 




Formalism was developed 
with the diagrammatical 
technique by Lifshitz 
and Pitaevskii 



the present work 




The connections between these ways have also been investigated and those 
have been confirmed in accordance together. 

Recently, the formalism of the nonequilibrium Green's functions and the 
nonequilibrium perturbative method have been applied widely to the vari- 
ous fields: condensed matter, plasmas, atoms and molecules, nuclear matter 
etc.. [5] Especially, the application to the problems of mesoscopic systems: 
the transports in quantum dots and quantum wire, has worked with great 
success[6-8]; for instance, electrical current and current noise are expressed 
in terms of the nonequilibrium Green's functions. For current, 

2eW r 

2(J.(z-l,z)) j dE[GU^,AE) - G<,_,AE)l (1.2) 

where W denotes the hopping matrix element. This reduces to the Landauer 
formula: [9] 

2(J) = ^ I dE[fL{E) - fR{E)]T{E), (1.3) 

where fi and are the Fermi distribution functions in the isolated left and 
right leads, respectively. For the current noise at zero-frequency from the 
autocorrelation function of the current [10]: 

2^°.(^ - 1,0 = -2 [-^j I dE[G<^,^{E) G>,^,^{E) 

-Gf__^ i_^^{E) Gfi„{E) 
-Gfm{E) Gf^^ i^^„{E) 
+GU,AE) GU,^{E) ]. (1.4) 

It reduces to the Khlus-Lesovik formula for one channel at zero-frequency for 
shot noise: [11] 

S^^ = 2(^-^^ jdE [ fn{E){l-ME)}r{E) 

+ h{E){l - h{E)}T{E) 

+ {fL{E)-fn{E)rT{E){l-T{E)}]. (1.5) 

Here, the transmission probability through a noninteracting system in Eqs. 
( 1.3 ) and ( 1.5 ) is written by 



T{E) = VL{E)G\^^{E)Vn{E)Gl,SE). 



(1.6) 



Tl and Tp> mean the coupling functions with the left and the right leads, 
respectively, and C and G"^ are retarded and advanced Green's functions, 
severally, as described in detail later. The expression for current noise also 
reduces to the Johnson-Nyquist noise for thermal noise: AQksT ( here, Q 
signifies conductance, and ks is the Boltzmann's constant and T denotes 
temperature ).[12] These are excellently compatible with the values observed 
by experiments. 

Nonetheless, the formalism of the nonequilibrium perturbation theory has 
not been completed yet; the perturbative methods with the diagrammatic 
technique remain to be clarified well. In particular, despite of that the defi- 
nition of the nonequilibrium Green's functions is given in time, the real-time 
perturbative expansion on the adiabatic theorem[13]— the basics of method 
of the perturbative expansion— is almost unknown. The general formalism 
to formulate arbitrary self-energy has still not been established and the con- 
nection with other formalism, for instance, the Matsubara imaginary-time 
perturbative formalism[14] is obscure. For this reason, the present work suc- 
ceeds to the Keldysh formalism and the diagrammatical method of Lifshitz 
and Pitaevskii to make progress in the formalism of the real-time perturba- 
tive expansion based on the adiabatic theorem. If this matrix form equation 
Eq. ( 1.1 ) works exactly as the Dyson's equation, this equation must be 
convertible into the Dyson's equations for retarded and advanced Green's 
functions. Accordingly, there must be the relations between self-energies. 
However, since the definition of the self-energies for this Dyson's equation in 
matrix form is not given, the self-energies drawn from perturbative expansion 
have not been made clear. 

In the present work, thus, the solution on the relations between self- 
energies is given. Using the solution, the retarded and advanced self-energies 
are derived from the self-energies in matrix form. Then, the derived retarded 
and advanced self-energies meet the conditions required as functions in time 
and the generally known relations on nonequilibrium Green's functions are 
fulfilled. Additionally, the formulated self-energies are in agreement with 
those derived by perturbative expansion in imaginary-time and analytical 
continuity. Thereby, the solution is in accordance with the generally known 
formation. As a result, it infers that the formalism of the nonequilibrium 
perturbation theory can be generalized. 



1.2 Kondo Effect 



The Kondo effect was discovered forty years ago [15]; the phenomenon of 
the minimum of the electrical resistivity in metals was explained in view of the 
interaction between conduction electrons and impurity by Kondo. After that, 
the Kondo physics was clarified from Landau's Fermi liquid theory [16], the 
renormalization group[17] and scaling[18]. Besides, the generalized Kondo 
problem with more than one channel or one impurity was proposed. [19] It 
has then been investigated in further detail. [20-23] Especially, the resistiv- 
ity has been expressed for the multichannel Kondo effect by the conformal 
field theoretical work, [20-22] in agreement with experiments, as mentioned 
in Chapter 4. 

Moreover, the Kondo effect in electron transport through a quantum dot 
was predicted theoretically at the end of 1980s[24-26]. After a decade, finally, 
this phenomenon was observed. [27] The Kondo effect was studied theoret- 
ically by use of the Anderson model. From scaling theory [18], the Kondo 
temperatures 

kBTK--De~^^"'^. (1.7) 

Here D means the band-width and the Coulomb interaction U , and F is the 
coupling function with leads, corresponding to the density of states for con- 
duction electron. The Kondo temperatures correspondent to the strength 
of the Kondo coupling decrease with increasing f/, as found from Eq. ( 1.7 
). The predictions from theoretical work using Anderson model were con- 
firmed experimentally. In the Kondo regime, the conductance was observed 
to reach the unitarity limit and the Kondo temperatures estimated from ob- 
servation [28] are in excellent agreement with the expression derived by the 
use of scaling theory for the asymmetric Anderson model [29]: 

where on-site energy is Eq. The perturbative approach, the Yamada-Yosida 
theory[30]— the perturbation theory for equilibrium based on the Fermi liquid 
theory[16] with the Matsubara imaginary-time perturbative method[14] is 
quite successful. 

Furthermore, the Kondo effect in a quantum dot was studied for nonequi- 
librium system where bias voltage is applied. [31] We have to know not only 



the Kondo effect but also the nonequihbrium state caused by bias voltage. 
The Yamada-Yosida theory was extended to nonequilibrium systems with the 
help of the Keldysh formalism and the Kondo effect in nonequilibrium system 
was studied. As the results, it was shown that for bias voltage higher than the 
Kondo temperatures, the Kondo resonance disappears in the spectral func- 
tion with the second-order self-energy of the Anderson model. [8] This results 
have been discussed little. After that, on experiments, it has been observed 
that the Kondo effect is suppressed when source-drain bias voltage is com- 
parable to or exceeds the Kondo temperatures. [32,33] The numerical results 
of the present work are also consistent with those. For the Kondo effect in 
nonequilibrium systems, it has been expected that the Kondo peak splits by 
bias voltage and that the two separated energy levels made in a quantum dot 
by the Kondo coupling act as two channels for two channel Kondo effect. In 
order to search that, every efforts have been done for many years. Such the 
phenomenon has however, never been observed for a quantum dot connected 
with two normal leads. In the present paper, it is discussed the reason why 
the Kondo peak is just broken and the Kondo peak splitting may not take 
place in simple two terminal systems with leads of the continuous energy 
states. 



Chapter 2 

Nonequilibrium Perturbation 
Theory 



A thermal average can be gained on the basis of the nonequihbrium per- 
turbation theory. [1,3,4,13,34-37] 

The perturbation theory is based on the adiabatic theorem ( called the 
Gell-Mann and Low's theorem[13] ). The Hamiltonian is given by 



Ti.0 and Tii are unperturbed and perturbed terms, respectively. Here, assum- 
ing that we can know only the state of Tio at t = — oo, thereby, initially 
Til = at t = — oo, so that the system is equilibrium and/or noninteracting 
state. The perturbation is turned on at t = — oo and introduced adiabat- 
ically. Then the perturbation is brought wholly into the system at t = 0; 
around t = 0, the system is regarded as stationary nonequilibrium and/or 
interacting state. After that, the perturbation is taken away adiabatically 
and disappears at t = oo. When the time evolution of the state is reversible, 
the state at t = oo can be expressed using the state at t = — oo by adding 
the phase factor. 

Now, let us consider the nonequilibrium state. If the time evolution of 
the state is irreversible for the nonequilibrium state, then, the state at t = oo 
is not well-defined; when the perturbation is removed entirely at t = oo, the 
state does not come back to the same state as at t = — oo. In this case, the 
ordinary perturbative method should be improved; the time evolution should 
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return to the well-defined state at t = — oo. Accordingly, the time evolution 
is performed along the real-time contour which starts and ends at t = —oo 
as illustrated in Fig. 2.1. It is the extension of the Gell-Mann and Low's 
theorem. [13] 



+branch 
-OO ^ 



3 



oo 

-oo ^ 

-branch 



Figure 2.1: time contour which starts and ends at t = — oo through t = oo 



2.1 S-matrix (S-operator) 

The Hamiltonian is given by Eq. ( 2.1 ). The time evolution in the 
interaction representation is expressed in terms of S-matrix by 

iP{t) = S{t,to)tP{to). (2.2) 

S-matrix S{t, to) is defined by 

Sit, to) = e^Wot/?ig-iH(t-to)/?ig-^Woto/?i^ (2.3) 

and has the following properties: 

S{to,to) = 1, (2.4) 
S{t,to)^ = S{to,t), (2.5) 

Sit,,t2)S{t2,ts) = 5(tl,t3), (2.6) 

^^aS|to) ^ ^^s{t,to). (2.7) 



Equation ( 2.7 ) can be solved formally by 



oo 1 

1 + E 



n=l 



" rt 



to 



dt„T 



to 



T 

S{to,t) = T 



exp i y dt Hi{t 



exp s — / (it 7Yi(t 



(2.8) 
(2.9) 



Here, the time ordering operator T arranges in chronological order and T is 
the anti time ordering operator which arranges in the reverse of chronological 
order. For the time evolution along the time contour as in Fig. 2.1, S- 
matrices, Eqs. (2.8) and (2.9) are required for paths on the (— ) branch 
from t = — oo to t = cxD and on the (+) branch from t = oo to t = —oo, 
respectively. 

In the same way, S-matrix for imaginary-time perturbative formalism is 
defined by 

S^{r, To) = e^«"/^-^("-"«)/'*e-^°"°/^ (2.10) 
and is also written by 

5.(r,ro) = l + E^hr / ^^i---/ dT„Tr niin)...niiT„) 

^-^nl \ h J Jto Jto L J 

expj^ £dT'ni{T)^y (2.11) 

These are the requisites for the Matsubara imaginary-time perturbative for- 
malism. [14] 



2.2 Matsubara Imaginary- Time Perturbative 
Formalism 



For thermal equilibrium, the statistical operator ( density matrix ) is 
written in Gibbs form for the grand canonical ensemble by 



e 



Qg = 



Tre-/3(w-M^) 



(2.12) 



13 = l/ksT. By rearranging Eq. ( 2.10 ), we have 



(2.13) 



By substitution of Eq. ( 2.13 ) into Eq. ( 2.12 ), the thermal average for 
equilibrium is obtained by 



By insertion of Eq. ( 2.11 ) in Eq. ( 2.14 ), the perturbative expansion is 
executed for the Matsubara imaginary-time perturbative formalism [14] us- 
ing the Bloch and De Dominicis's theorem[38-41]. The Matsubara Green's 
function is defined by 



The functions in terms of the Matsubara Green's function in imaginary- 
time are converted into those in the Matsubara frequency by the Fourier 
transformation for the Matsubara Green's function: 



The Matsubara frequency, Un= (2n + l)7r//3 for fermion and ujn= 2mT / (5 for 
boson, (n = 0, ±1, ±2, ±3 ... ); thereby, the Matsubara Green's function is 
periodic. 

After that, the analytical continuity is performed. For fermion, when the 
Taylor expansion for the function e^^ + 1 around poles z = on imaginary 
axis is done, then, e^^ +1 = 0, hence, z^ = (2n + [3, and as approx- 
imation, e^^ + lRi/3e^^ {z — z^). The residue theorem yields the conversion 
of sum of the functions in the Matsubara frequency into the contour integral 



(...) 



Tr[e-^(^-^^)...] _ Tr[e-^(^°~^^)^^(/3/i,0)...] 
Tr[e-/3(w-A'^)] ~ Tr[e-/5(^o-/.^)^^(/5^,o)]. 



(2.14) 



^(r)^- (T,d(r)c/t(o)). 



(2.15) 




(2.16) 



by 




(2.17) 



It should be noted that the contour C in integral surrounds the poles on 
imaginary axis. After that, for example, when g{z) = l/{z — e^), then, the 
contour integral with contour C enclosing e^, a pole on real axis is executed 
by 

where /(e^) is the Fermi distribution function. 

In the same way, in boson case, for poles z = z^ on imaginary axis, 
e^^ —1 = 0, i.e. = 2n7ri/j3, and approximately, e^^ — l^(3e^^ {z — z^), 
we have 

^E^?M = -— /-f^ (2.19) 

Then, by changing the contour integral around poles on imaginary axis 
into the contour integral parallel with real axis ( e.g. E±i6 ), we have 
the functions written in terms of retarded and advanced Green's functions 
in energy. For high-order perturbation theory, the analytical continuation 
is so complicated. As the method of analytical continuation, Eliashberg's 
method[42] is known. [4,41] 



2.3 Nonequilibrium Perturbative Formalism 

2.3.1 Nonequilibrium Real-Time Perturbative Formal- 
ism 



For nonequilibrium. Equation ( 2.12 ) is not exact. We should note that 
there are no specific limitations upon the statistical operator. The von Neu- 
mann's statistical operator is expressed independently of whether the states 
are at thermal equilibrium or nonequilibrium by[35,37] 

Qs{t) = T.\^s{t))Pn.{ms{t)\ (2.20) 

m 

in the Schrodinger representation. Here, Pm is probability that the system 
is in state m and \ms{t)) is the state in the Schrodinger representation. It 



satisfies the Liouville equation by 

zh^ = [n,gs]. (2.21) 
The statistical operator in the interaction representation is given by 

g(t) = e^^»*/^^?s(^)e-^^»*/^ (2.22) 
and also obeys the Liouville equation by 

= [ni,g]. (2.23) 

As a matter of course, 

te(0) = Qh{0) = (2.24) 

where QHit) is statistical operator in the Heisenberg representation. The 
time evolution is written by means of S-matrix by 

g{t) = S{t,to)Q{to)S{to,t). (2.25) 

These are the properties of the von Neumann's statistical operator. Although 
the von Neumann's statistical operator Eq. ( 2.20 ) is not an explicit expres- 
sion, it is still considered to be the same type as Eq. ( 2.14 ) so as to execute 
the perturbative expansion. 



The thermal average for nonequilibrium is drawn in view of the analogy 
with the imaginary-time perturbative method for the Matsubara Green's 
function using S-matrix and the time evolution of the statistical operator, 
Eq.( 2.25 ). Thus the thermal average of the operators in the Heisenberg 
representation at t = can be brought, for example by[3,4,34,37] 

{TAit)B{t)) 

= Tr 



Tr 



gH{0)TA{t)B{t 
^(-oo)5(-oo, 0)TA(t)5(t')5(0, -oo] 



Tr |^f)(— oo)5(— oo, oo) {T5(oo,-oo)i(r)5(t'-)} 

j:j:^^i-n} (t) Lj'^-Lj'-L 

{ {t w,(«f ). . .Hi(C)} {th,(«;-). . .n,(i-)A{t-)B{t'-)} ) 



n=l m=l 
X 



(2.26) 



where (. . .)^^ means Tr[^(— oo). . .]. A denotes an arbitrary operator in the 
interaction representation. Here, A(t)= S{0,t)A{t)S(t,0) and ^_h'(0)= ^(0) 
= S{0, —oo)g{—oo)S{—oo,0). In this case, the operators A{t~) and B{t ~) 
in Eq.( 2.26 ) both are on the (— ) branch from t = — oo to t = oo of the time 
contour in Fig. 2.1. 

Using the last expression in Eq.( 2.26 ), the real-time perturbative expan- 
sion is executed diagrammatically by the help of the Wick's theorem. On the 
diagrammatical perturbative expansion, the diagrammatical methods writ- 
ten by Lifshitz and Pitaevskii. [4] are extended. 



+ 



Figure 2.2: The extension of the diagrammatic technique of Lifshitz and 
Pitaevskii. S^^^ brought as the sum over terms going from time (— ) to 
time (— ) by way of times (— ) and (+). 
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Figure 2.3: S^"*-^^ brought as the sum of terms passing by way of times, 
(-)(-), (+)(+), (-)(+) and (+)(-). 

There the summation over terms in all times is taken. For example, the 
third-order self-energy E^^^ is brought as the sum over terms going from 
time (— ) to time (— ) by way of times (— ) and (+) as illustrated in Fig. 2.2. 
See Fig. 2.1 again. It should be noted that the times (— ) are on the (— ) 
branch from t = —oo to t = oo of the time contour in Fig. 2.1 and the times 



(+) are on the (+) branch from t = oo to t = — oo of the time contour. In 
addition, for the fourth-order self-energy, Ti^~^^\ the sum of terms passing 
through times, (—)(—), (+)(+), (~)(+) and (+)(—) is taken, as shown in 
Fig. 2.3. 

2.3.2 Schwinger-Keldysh Formalism 



The Dyson's equation for the Keldysh formahsm is given by 



G 



S G, 



where 



G 



G- 

G> 



G< 

G+- 





" s— 


s< 




s> 









The nonequihbrium Green's functions are defined by 

G-{t^M)= -i{Td{h)d^{t2)), 

G>{hM)= -i{d{h)d\h)), 
G<{UM)=i{d\h)d{t{)). 



(2.27) 



(2.28) 
(2.29) 
(2.30) 
(2.31) 



As shown in Eq. ( 2.26 ), the operators in Eq. ( 2.28 ) both are on the (— ) 
branch from t = —oo to t = oo of the time contour in Fig. 2.1. In the same 
way, those in Eq. ( 2.29 ) both are on the (+) branch from t = ootot = — oo 
of the time contour, and those in Eq. ( 2.30 ) are on the (+) branch and on 
the (— ) branch, respectively. In Eq. ( 2.31 ), they are on the (— ) branch 
and on the (+) branch, respectively. 



Additionally, retarded and advanced Green's functions are defined by 

G'it,,t2)= -teih-t2){{d{h)J\h)}), (2.32) 
G''ih,t2)=ieih - h){{dit,)J\h)}). (2.33) 

Here, the curly brackets signifies anticommutator. The Dyson's equations 
for retarded and advanced Green's functions are given by 

G' = / + / G', (2.34) 
G" = + G". (2.35) 



As the necessity to Eqs. ( 2.34 ) and ( 2.35 ), the self-energies S'' and 
must be retarded and advanced functions in time, respectively. 



In accordance with the ordinary procedure of nonequilibrium perturbative 
formalism, [3,4, 34,35,37] for the Dyson's equation for the Keldysh formalism, 
Eq. ( 2.27 ), the transformation is carried out: by use of 

1 -1 ' 
1 1 



1 

71 



G 



and 



G> 



G< 



LGL^ 







G" 







From the definition of Green's functions, 



in other words, 



G' 



G^{t) 
G\t) 



G- 
G- 



G< 
G> 



G> -G^ 
G< -G^ 



[G>{t)-G<{t)]e{t), 

[G<{t)-G>{tm-t). 



G is called the Keldysh Green's function. For the self-energies part, the 
following relationship is required: 



VL{t) 



s— (t) + S<(t) = 
(t) + S>(t) = 

s-(t) + s++(t) 



-S++(t)-S>(t), 
-S++(t)-S<(t), 
-S<(t)-S>(t). 



(2.36) 
(2.37) 
(2.38) 



The above is known in general. Here, it is uncertain whether or not the 
requirements of self-energies as functions in time are fulfilled. It is because 
the definition of the self-energies in matrix form is not given, as mentioned 
earlier. 



The solution is given from the present work. [43] That is explained as 
follows: in the right-hand sides of Eqs. ( 2.36 ), ( 2.37 ) and ( 2.38 ), the 
functions, S (t), S+"''(t) and are deduced from perturbative 

expansion with the Wick's theorem. At this time, the directions in time are 
not defined in the functions S (t), S'^(t), T,~^~^{t) and T,^{t). This is due 
to the following functions in dependence upon time: 

g--{t) = e{t)g>{t) + e{-t)g<{t), (2.39) 
g++it)=eit)g<it) + e{-t)g>it). (2.40) 

For this reason, for the right-hand sides of Eqs. ( 2.36 ), ( 2.37 ) and ( 2.38 
), the directions in time must necessarily be taken into consideration. 

As the main point, the terms in the right-hand sides are taken for sum of 
retarded and advanced terms: 



S^(t) = [S"-(t) + S<(t)]^(t) + [S-"(t) + S<(t)]e(-t) 

= -[S++(t) + S>(t)]^(t)-[S++(t) + S>(t)M-t), (2.41) 



E^{t) = [s--(t) + s>(t)]e(t) + [s-(t) + s>(t)]^(-t) 

= -[S++(t) + S<(t)]^(t) - [S++(t) + S<(t)]^(-t), (2.42) 



n{t) = [s--(t) + s++(t)]^(t) + [s--(t) + s++(t)]^(-t) 

= -[S<(t) + S>(t)]e(t)-[S<(t) + S>(t)]e(-t). (2.43) 



For self-energy functions derived by the present perturbative expansion 
via the Wick's theorem, the following relations are found: 

s— (t)^(t) = -j:>{t)e{t), (2.44) 

S++(t)^(t) = -S<(t)^(t), (2.45) 
S-"(t)^(-t) = -S<(t)^(-t), (2.46) 
S++(t)e(-t) = -S>(t)^(-t); (2.47) 



these relations have never been known in generaL When these are substituted 
into Eqs. ( 2.41 ), ( 2.42 ) and ( 2.43 ), then, the advanced term of Eq. ( 
2.41 ) and the retarded term of Eq. ( 2.42 ) are canceled: 

[S~(t) + S<(t)]^(-t) = -[S++(t) + S>(t)]^(-t) = 0, (2.48) 
[S-(t) + S>(t)]e(t) = + S<(t)]^(t) = 0, (2.49) 

so that Equations ( 2.41 ) and ( 2.42 ) reduce to 

E'it) = [S--(t) + S<(t)]e(t) = -[S++(t) + S>(t)]e(t), (2.50) 
S'^(t) = [S--(t) + S>(t)]^(-t) = -[S++(t) + S<(t)]^(-t). (2.51) 

The retarded and advanced self-energies are acquired as retarded and ad- 
vanced functions in time, respectively; they are the requirements. In addi- 
tion. Equation ( 2.43 ) is certainly reproduced. 

Then, it leads to 

Z^{t) = [S<(t)-S>(t)]e(t), (2.52) 
S"(t) = [S>(t)-S<(t)]^(-t); (2.53) 

they consist with the expressions in the review by Rammer and Smith. [44] 
The formalism in the review of Rammer and Smith is different from the 
nonequilibrium perturbative expansion, so that the present method is con- 
firmed to connect with the other formalism. By performing the Fourier trans- 
formation for Eqs. ( 2.52 ) and ( 2.53 ), we have 

S'^(E) - S'^(E) = S<(E) - S>(E); (2.54) 

it proves that the relation stands. Since Equation ( 2.54 ) is widely known. 
Equations ( 2.52 ) and ( 2.53 ) should generally hold as functions in time. 
Besides, 

G< = (1+G'^S'^)^<(1 + G"^S'^)-G'^S<G^ (2.55) 
G> = (1+G"^S'^)^>(1 + G"^S'^)-G'^S>G'^ (2.56) 

and 

G^ = (1 +G"S")^-^(1 + G'^S'') (2.57) 

still work. 

As mentioned above, the present solution is in accordance with the gen- 
erally known relations. It indicates that the present solution has validity. 



Chapter 3 

Expressions of Self-Energy for 
Anderson model 



3.1 Anderson model 



We consider equilibrium and nonequilibrium stationary states. Nonequi- 
librium state is caused by finite bias voltage, that is, the difference of chemical 
potentials; after bias voltage was turned on, long time has passed enough to 
reach stationary states. Since the states are stationary, the Hamiltonian has 
no time dependence. The system is described by the Anderson model linking 
to leads. The impurity ( the quantum dot ) with on-site energy Eq and the 
Coulomb interaction U is connected to the left and right leads by the mixing 
matrix elements, vl and vr. The system is illustrated below. 



U=0 



Ut^O 



Left Lead 



Quantum 
Dot 



U=0 



Right Lead 



The Anderson Hamiltonian is given by 
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- vlAcl. + H.c.) - vnidicR^ + H.c). (3.1) 

d) (d) is creation (annihilation) operator for electron on the impurity, and 
and c|j (c^ and cr) are creation (annihilation) operators in the left and 
right leads, respectively, a is index for spin. The chemical potentials in the 
isolated left and right leads are fiL and fiR, respectively. The applied voltage 
is, therefore defined hj eV = fiL — fin- 

We consider that the band- width of left and right leads is large infinitely, 
so that the coupling functions, Tl and Tr can be taken to be independent 
of energy, E. On-site energy Eq is set being canceled with the Hartree term, 
i.e. the first-order contribution to self-energy for electron correlation, as 
mentioned later. 

Accordingly, the Fourier components of the noninteracting ( unperturbed 
) Green's functions reduce to 

fAE) = ^ (3.2) 

,'(E) = ^ (3.3) 

where F = (Fj^ + F/j)/2. Hence, the inverse Fourier components can be 
written by 

/(t) = -t9it)e-^\ (3.4) 

g'^it) = i9{-t)e^\ (3.5) 

In addition, by solving the Dyson's equation Eq. ( 2.27 ), we have 

g<{E) = g^\E) [z/i(i?)Fi + tfR{E)TR]g\E), (3.6) 

g>{E) = g^\E) WdE) - 1)Vl + iUr{E) - l)TR]g'^{E). (3.7) 

fi and Jr are the Fermi distribution functions in the isolated left and right 
leads, respectively. By Eqs. ( 3.6 ) and ( 3.7 ), the nonequilibrium state is 
introduced effectively as the superposition of the left and right leads. In this 
case, the effective Fermi distribution function can be expressed by[8] 

f (Tp\ fL{E)TL + fR{E)TR 

JcS[E) = — — . (3.8j 

i L + i i? 

This effective Fermi distribution function is reasonable because it is con- 
sidered that leads have the continuous energy states and the two chemical 
potentials are not two localized states. 



3.2 Self-Energy 



The expressions for self-energies of the Coulomb interaction of the Ander- 
son model are formulated by the method of the nonequilibrium perturbation 
theory based on the adiabatic theorem, explained in Chapter 2. Practically, 
the perturbative expansion is done with respect to the Coulomb interaction 
term of Eq. ( 3.1 ) using the last expression of Eq. ( 2.26 ) in view of the 
Dyson's equation Eq. ( 2.27 ) and each diagram. In S (t), 
and S>(t) obtained by the present solution. Equations (2.44)-(2.47) are sat- 
isfied. Every formulated Ti^'{E) and are in the relation of the complex 
conjugate each other. 



3.2.1 First-Order Contribution 



For the Anderson model, the first-order contribution, the Hartree term 
can be written by 

E^WfE) = = U{n) = U f —G<(E), (3.9) 

J 27ii 

where (n) is charge density. 



3.2.2 Second-Order Contribution 



Jo 



The second-order self-energies are expressed by 

g>{t,)g>{t,)g<{-t^) 

-9^{ti)g<{ti)g>{-h] 

9Hti)9^ih)g<{-t,) 
+g<{t,)g>{t,)g^{-t,) 



dtie 



lEti 



J —oo 



iEti 



g<{h)g<{U)g>{-U) 
-9^{h)g>{h)g<{-h] 

9^{ti)9^ih)g<i-t,) 

+9Hti)9Hti)9^{-ti) 
+g<{t,)g>{t,)g^{-t,) 



(3.10) 



(3.11) 



Here g^{t) = g''(t) + g^it), that is, g'^it) = g''{t) = -ie{t)e-^^ for t>0 and 
g-[t) = g^it) = ie{-t)e^^ for t < 0. Additionally, g<{t) and g>{t) are the 
inverse Fourier components of Eqs. ( 3.6 ) and ( 3.7 ). Figure 3.1 shows 
the diagram for the second-order self-energy. These expressions for equilib- 
rium agree exactly with those deduced from the Matsubara imaginary-time 
perturbative expansion for equilibrium and analytical continuity by Zlatic et 
a/. [45]. As shown numerically later, the second-order contribution coincide 
with those brought out by Hershfield et a/. [8]. 

In the symmetric equilibrium case, the asymptotic behavior at low energy 
is expressed by 

the exact results based on the Bethe ansatz method. [46,47] 

► 



Figure 3.1: The second-order self-energy 



3.2.3 Third-Order and Fourth-Order Contributions 



There are two kinds of the third-order contributions as illustrated in Fig. 



3.2. 
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g<{-h)g>ih-t2)g>ih-t2) 
-9^{-ti)g''{ti-t2)g<{h-t2) 



X 



g^{t2)g>{t2) + g<{t2)gHt2) 



(3.13) 



/O poo 
-oo J —oo 



g>{-U)g<{U-t2)g<{U-t2) 
-9'^{-ti)g>{ti-t2)g>{ti-t2) 



g^{t2)g>{t2)+g<{t2)gHt2] 



(3.14) 



-r{3) 



{E) 



X 



dti / dt2e 



iEti 



g>{U)g>{U-t2)g^{t2-h) 
-g<{ti)g<{h-t2)g>{t2-h] 



gHt2)g<{-t2)+g<it2)gH-t2) 



(3.15) 



y^a(3) 



dti 



due' 



g<{ti)g<{ti-t2)g>{t2-h) 
-g>{ti)g>{h-t2)g<{t2-h] 



gHt2)g<{-t2) + g<{t2)g^ht2) 



(3.16) 



Equations (3.13)-(3.16) for equilibrium state agree exactly with those de- 
rived from the Matsubara imaginary-time perturbative expansion for equi- 
librium and analytical continuity by Zlatic et a/. [45]. As mentioned later, 
it is numerically confirmed that the third-order contribution is canceled for 
the symmetric Anderson model; this is compatible with both the results de- 
duced from the Yamada-Yosida theory[30, 47,48] and gained on the basis of 
the Bethe ansatz method[46]. 





Figure 3.2: Two kinds of the third-order self-energies: E^p*(Left) and S^^^ 



Furthermore, the fourth-order contribution to the self-energy is formu- 
lated. ( See Appendix. ) The twelve terms for the proper fourth-order 
self-energy can be divided into four groups, each of which comprises three 
terms, corresponding diagrams Fig. 3.3 (a)-(c), (d)-(f), (g)-(i) and (j)-(l), 
respectively. 
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Figure 3.3: The twelve terms for the proper fourth-order self-energy divided 
into four groups: (a)-(c), (d)-(f), (g)-(i), and (j)-(l). 



For symmetric Anderson model at equilibrium, the asymptotic behavior 
at low energy is approximately in agreement with those based on the Bethe 
ansatz method[46]: 

(3.17) 

The fourth-order contribution has not been clarified well. In particular, 
the behavior for nonequilibrium state has almost been unknown. In Chapter 
4, the numerical results are shown and discussed. 



Chapter 4 



Numerical Results and 
Discussion 



4.1 Self-Energy 



The third-order terms, Eqs. ( 3.13 )-( 3.16 ) are canceled under electron- 
hole symmetry not only at equilibrium but also at nonequilibrium: 



As a consequence, the third-order contribution to self-energy vanishes in the 
symmetric case. It is consistent with the results of Refs. [30, 47, 48] based 
on the Yamada-Yosida theory that all odd-order contributions except the 
Hartree term become null for equilibrium in the symmetric single-impurity 
Anderson model; probably, it is just the same with nonequilibrium state. 
On the other hand, the third-order terms contribute to the asymmetric sys- 
tem where electron-hole symmetry breaks and furthermore, the third-order 
terms for spin-up and for spin-down contribute respectively when the spin 
degeneracy is lifted for example, by magnetic field. For the fourth-order con- 
tribution, three terms which constitute each of four groups contribute equiv- 
alently under electron-hole symmetry. Moreover, to the asymmetric system, 
the terms brought by the diagrams of Fig. 3.3(a) and (b) contribute equiva- 
lently and the terms by the diagrams of Fig. 3.3 (j) and (k) make equivalent 
contribution, and the rest, the eight terms contribute respectively. Further, 
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Figure 4.1: The second-order self-energy for the symmetric Anderson model 
at U /T = 1.0 and zero temperature. Left:real part and Right: imaginary part 
at equilibrium ( solid line ), eV/T = 1.0 ( thin sohd line ) and eV/T = 2.0 ( 
dashed line ). 

the twenty-four terms for spin-up and spin-down take effect severally in the 
presence of magnetic field. 

The second-order and the fourth-order contributions to self-energy for 
zero temperature symmetric Anderson model are plotted in Fig. 4.1 and 
in Fig. 4.2, respectively. Equation ( 3.12 ) represents the curves around 
E = denoted by solid line in Fig. 4.1, and Equation ( 3.17 ) represents 
approximately those shown in Fig. 4.2. The curves of the second-order self- 
energy shown in Fig. 4.1 are identical with those of expressions derived by 
Hershfield et al.[8]. In comparison of Fig. 4.2 with Fig. 4.1, it is found that 
the fourth-order contribution for equilibrium has the same but narrow curves 
at low energy with those of the second-order contribution. In addition, the 
broad curves are attached at high energy for the fourth-order self-energy. ( 
The higher contribution is, the more should the curves oscillate as a function 
of energy. ) When the voltage, eV/T exceeds ~2.0, the behavior of curves of 
self-energy changes distinctly and comes to present a striking contrast to that 
for the second-order contribution. Especially, the curve for the imaginary 
part of the fourth-order contribution rises up with maximum a.t E = 0. 
On the other hand, for the second-order contribution, a valley appears with 
minimum at energy of zero— it is quite the contrary. Moreover, from these 
results, it is expected that the sixth-order contribution to imaginary part of 
self-energy has minimum a.t E = 0. By reason of these, the perturbative 
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Figure 4.2: The fourth-order self-energy for the symmetric Anderson model 
aXU/T = 1.0 and zero temperature. Left:real part and Right: imaginary part 
at equilibrium ( solid line ), eV/T = 1.0 ( thin solid line ), and eV/T = 2.0 
( dashed line ). The fourth-order contribution for equilibrium has the same 
but narrow curves at low energy with those of the second-order contribution. 

expansion is hard to converge for eV/T > ~2.0, as mentioned later. 

4.2 Current Conservation 

In this section, the problem on the current conservation is described be- 
low. In Ref . [8] , it is shown that the continuity of current entering and leaving 
the impurity stands exactly at any strength of U within the approximation 
up to the second-order for the symmetric single- impurity Anderson model. In 
comparison of Fig. 4.2 with Fig. 4.1, it is found that curves of fourth-order 
self-energy have the symmetry similar to those of the second-order. From 
this, it is anticipated that the current conservation are kept perfectly with 
approximation up to the fourth-order in the single-impurity system where 
electron-hole symmetry holds. The continuity of current can be maintained 
perfectly in single-impurity system as far as electron-hole symmetry stands. 
On the other hand, current comes to fail to be conserved with increasing U 
in asymmetric single-impurity case and in two-impurity case. 



4.3 Spectral Function 

4.3.1 For Second-Order Self-Energy 
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Figure 4.3: The spectral function with the second-order self-energy at f//r = 
10.0 for the symmetric Anderson model at equilibrium ( solid line ), eV/T = 
1.0 ( thin solid line ) and eV/T = 2.0 ( dashed line ). 

The spectral function with the second-order self-energy is widely known. 
[8] [47] It is plotted for U/T = 10.0 and zero temperature in Fig. 4.3. For 
equilibrium, the Kondo peak at energy of zero is very sharp and the two-side 
broad peaks appear at ±U /2. The curve for eV = is identical with that 
shown in Ref. [47]. As eV becomes higher than the Kondo temperatures, 
A;^Tx[49], the Kondo peak becomes lower and is lost finally, while the two- 
side broad peaks rise at -E^ ±U /2.[8] 

4.3.2 For Self-Energy up to Fourth-Order for Equilib- 
rium 

Figure 4.4 shows the spectral function with the self-energy up to the 
fourth-order for equilibrium and zero temperature. With strengthening U, 
two-side narrow peaks come to occur in the vicinity oi E = ±U/2 in addition 





E/r 

Figure 4.4: The spectral function with self-energy up to the fourth-order at 
equilibrium for the symmetric Anderson model at f//r = 3.5 ( dashed line ) 
and U/T = 5.0 ( solid line ). 



to the Kondo peak. At U large enough, the Kondo peak becomes very acute 
and two-side narrow peaks rise higher and sharpen; the energy levels for the 
atomic limit are produced distinctly. As before, the fourth-order self-energy 
has the same but narrow curves as functions of energy with those of the 
second-order and those curves make the peaks a.t E = ±U/2. 

For the present approximation up to the fourth-order, the Kondo peak 
at £^ = reaches the unitarity limit and the charge, (n) corresponds to 1/2, 
that is, the Friedel sum rule is correctly satisfied [50]: 

p{Ef) = sin2(7r(n))/7rr. (4.1) 

where p{Ef) is the local density of states at the Fermi energy. 

Here, the discussions should be made on the ranges of U in which the 
present approximation up to the fourth-order stands. From the results, it is 
found that the approximation within the fourth-order holds up to U/T ~5.0 
and is beyond the validity for U /T> ~6.0. In such the higher-order 

terms are required. 



4.3.3 For Self-Energy up to Fourth-Order for Nonequi- 
librium 
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Figure 4.5: The spectral functions with self-energy up to the fourth-order 
at eV/T = 0.5 ( Left ) and eV/T = 1.0 ( Right ) for the symmetric Anderson 
model. U/r = 3.5 ( dashed line ) and U/T = 5.0 ( solid line ). 



Next, the results for nonequilibrium and zero temperature are shown. 
The expression for the Friedel sum rule, Eq. ( 4.1 ) does not stand for 
nonequilibrium, since the charge cannot be expressed with respect to the 
local density of states. All the same, the Kondo peak reaches the unitarity 
limit and (n) = 1/2 in the symmetric and noninteracting case. 

The spectral functions with the self-energy up to the fourth-order are 
plotted for eV/T = 0.5 and eV/T = 1.0 in Fig. 4.5, respectively. The 
Kondo temperatures lessen with the rise of U, as in general known. Approx- 
imately [49], 



kBTK-^^e-^/'^^-"^/'^, (4.2) 

As the estimation from Eq. ( 4.2 ), ksTx/T ~0.5 for U/T = 3.5 and kBTx/T 
~0.3 for U/r= 5.0. When U is strengthened and eV exceeds ksTK, the 
Kondo peak for eV/T = 0.5 falls in and instead, the two-side narrow peaks 
remain to sharpen in the vicinity of E = ±U/2. For eV/T = 1.0, the Kondo 
peak becomes broad and disappears for U large enough. The two-side peaks 



is generated small in the vicinity oi E = ±U/2. The Kondo resonance is quite 
broken for bias voltage exceeding the Kondo temperatures; this accords with 
the experimental results of two terminal systems that the Kondo effect is 
suppressed when source-drain bias voltage is comparable to or exceeds the 
Kondo temperatures, eVyksTx -[32,33] 

For eV/r > ~2.0, the Kondo peak does not lower even when eV is much 
larger than kBTx- The perturbative expansion is hard to converge on account 
of the imaginary part of the self-energy for eV/T > ~2.0, as described before; 
thereby, the higher-order contribution is probably needed for high voltage. 
Additionally, for high bias voltage, the picture of the nonequilibrium state 
represented as the superposition of the two leads is hard to stand. 

We have to consider not only the Kondo effect but also the nonequi- 
librium state caused by bias voltage. In the present work, nonequilibrium 
state is represented as the superposition of the two leads. As is known, this 
works well as expression of nonequilibrium state. In this connection, the 
Khlus-Lesovik formula [11], the expression for current noise is drawn from 
the present picture [10] and is quantitatively consistent with experiments in 
ballistic systems. Thereby, the present picture can be valid as description of 
nonequilibrium state. The effective Fermi distribution function, Eq. ( 3.8 
) is found qualitatively similar to that for finite temperatures fxiE). For 
example, for finite temperatures Tt^O, 

it is indicated that there are thermal charge fluctuations. In case of the 
effective Fermi distribution function, 

/eff(^)[l-/efr(i?)]^0 

for finite voltage eV^Q. From this analogy in the Fermi distribution function, 
it is inferred that there are nonequilibrium fluctuations similar to thermal 
fluctuations. Because of the effective Fermi distribution function, not only 
for the second-order but also for the fourth-order, the Kondo resonance is 
destroyed, qualitatively the same as for finite temperatures, [51] as is observed 
in experiments of two terminal systems. This is because the energy states in 
the leads are not discrete but continuous. The present results indicate that 
the Kondo resonance splitting due to bias voltage does not take place for 
simple two terminal systems. 



Over a decade, it has been expected that the Kondo resonance sphts 
by bias voltage and work as two channels for two channel Kondo effect, so 
that many endeavours have been made to seek the Kondo peak splitting. 
Recently, it has been reported that the two channel Kondo effect is observed 
in quantum dots system. [52] On the experiment, it is not observed for two 
infinite leads system which has continuous energy state. That occurs for the 
system that a large dot with closely discrete energy levels is added to two 
infinite leads system where a large dot works as finite size reservoir and two 
infinite leads system acts as an infinite reservoir. The observed differential 
conductance depends on bias voltage and temperatures; this characteristic 
shows the functions of bias voltage and temperatures for two channel Kondo 
effect, derived by the theoretical work with conformal field theory. [20-22] 



Chapter 5 
Summary 



In the present work, the solution is proposed so as to make progress in the 
Schwinger-Keldysh formahsm of nonequihbrium perturbation theory. Using 
the solution, the expressions for self-energies are derived from the pertur- 
bative expansion in real-time through the Wick's theorem. On derivation, 
the expressions drawn via perturbative expansions can be taken for sum of 
retarded and advanced terms, and it is demonstrated that the advanced term 
in the retarded self-energy and the retarded term in the advanced self-energy 
vanish. As the consequence, the retarded self-energy is obtained as retarded 
function in time and the advanced self-energy is expressed as advanced func- 
tion. These expressions for self-energies agree with those acquired by execut- 
ing perturbative expansion in imaginary-time and analytical continuity; it is 
proven that the method of nonequihbrium perturbation theory can be linked 
with that of perturbation theory of thermal Green's functions. It is indicated 
that the Schwinger-Keldysh formalism works appropriately and that the for- 
malism can be united with the other valid methods. The Schwinger-Keldysh 
formalism should be clarified and generalized furthermore. 

As the numerical results, the Kondo peak disappears as bias voltage ex- 
ceeding the Kondo temperatures. Because of the analogy of the effective 
Fermi distribution function for nonequihbrium with that for finite tempera- 
tures, the present result is qualitatively similar to that for finite temperatures. 
This characteristic appears in the experiments of two terminal systems. Con- 
sequently, the split of the Kondo peak by finite bias voltage may not occur 
in simple two terminal systems by reason of that the leads have continuous 
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energy state. 

The nonequilibrium systems in the presence of bias vohage are too com- 
phcated and of difficult access. Apart from the nonequilibrum perturbation 
theory as the present work, the statistical physics approach to the nonequi- 
librium systems is well-known as nonequilibrium statistical mechanics and 
nonequilibrium statistical thermodynamics by Zubarev; nonequilibrium sta- 
tistical operator has also been studied as the generalized Gibbs operator ( 
Gibbs operator is given by Eq. ( 2.12 ) ).[53] To the problem on the nonequi- 
librium systems by bias voltage, various approaches of both the statistical 
physics and phenomenology are also required. As for the Kondo physics, the 
various problems about the Kondo effect in nonequilibrium systems and the 
multichannel Kondo effect remain to be clarified more. 



Chapter 6 
Appendix 



6.1 Appendix A 

Fourth-Order Contribution to Self-Energy 

The twelve terms for the fourth-order contribution can be divided into four 
groups, each of which is composed of three terms. The four groups are 
brought from diagrams denoted in Fig. 3.3 (a)-(c), (d)-(f), (g)-(i), and (j)- 
(1), respectively. The terms for the diagrams illustrated in Fig3.3 (a) and (b) 
are equivalent except for the spin indices and expressed by 
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Additionally, Figure 3.3(c) shows the diagram for the following terms: 
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Next, the terms brought from diagram in Fig. 3.3(d) are expressed by 
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The terms for diagram in Fig. 3.3(e) are written by 
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9^{h-t2)g>{h-t2)g<{h-U) 
-9^{ti-t2)g<{ti-t2)g>{h-h) 

9^it2-h)g>{-h)g<{h) 
-g^{h-h)g<{-h)g>{h) 



X g^it2)sgn{ts) 



(A.7) 



rO roc roc 

E^(4)(E) = f/M dti dt2 dt^e'^'^ 



X 



g''ih-t2)g<ih-t2)g>ih-h) 
-g''{ti-h)g>{h-h)g<{h-h) 

g''{h-h)g>{-ts)g<{ts) 
-g^it2-h)g''i-h)g>{h) 



X g^it2)sgn{ts) 



(Ai 



In addition, Figure 3.3(f) denotes the diagram for the following terms: 



dti / dt2 / dts e 



JEti 



g''{h-t3)g>{ti-t2)g<{t2-ti) 
-g^{ti~h)g<{ti-t2)g>{t2-h) 

g<{h)g<{h)g>{t2-h) 
-g^{t3)g^it3)g-^{t2-h) 



X ^^(-^2) sgn(t3) 



(A.9) 



Ef\E) = I dt, I dt2 I dhe'""'' 



X 



g^ih-t;)g<{t,-t2)g>{t2-h) 

-g^{h-h)g>{ti-t2)g<{t2-h) 

g<{h)g<{h)g>{t2-h) 
-g^{h)g>{h)g<{t2-h) 



X ^^(-t2)sgn(t3) 



(A.IO) 



Next, the terms formulated from diagram illustrated in Fig. 3.3(g) are ex- 
pressed by 



Ef)(E) 



X 



dt-2 



dU e 



iEU 



g>{ti)g>{h-t2-h)g<{t2-h) 

-g<{t,)g<{t,-t2-h)g>{t2-h) 

g''{t2 + t3)g>{h)g<{-h) 
-g^{t2 + h)g<{h)g>{-h) 



X ^^(-t2)sgn(t3) 



(A.ll) 



/O poo POO 

dti / dt2 / dh e*^*i 
-OO J —OO J —CO 



X ^^(-t2)sgn(t3) 



-g>{t,)g>{t^~t2-h)g<{t2-h) 

g^it2 + h)g>{h)g<{-h) 
-g^{t2 + h)g<{h)g>{-h) 



Figure 3.3(h) illustrates the diagram for the following terms: 



T.'^'^\e) = I dti I dt2 I dhe 



AEti 



g<{h)g<{h-t2-h)g>{t2-h) 
-g>ih)g>{h-t2-h)g<{t2-h) 

g^ih)g>ih)g<i-t2-h) 
-g<{t;)g<{t;)g>{-t2-h) 



x^^(^2)sgn(t3) 



(A. 12) 



(A.13) 



Ef)(i?) 



X 



dti 



dU 



g>{ti)g>{ti-t2-h)g<{t2-U) 
-g<{t,)g<{t,-t2-h)g>{t2-h) 

g>{h)g>{h)g<{-t2-ts) 
-g''ih)g<{h)g>i-t2-h) 



x^^(t2)sgn(t3) 



(A. 14) 



Besides, the terms formulated from the diagram in Fig. 3.3(i) are written by 

/ A\ P'^ P'^ 

Yr^^\E) = / dti / dt2 / rftae'^*^ 

JO J —OO J —OO 



X 



g>{-ti)g<{h-t2-h)g<{ti-t2) 
-g^i-ti)g>ih-t2-ts)g>ih-t2) 

g''it2 + h)g>{h)g<{-h) 
-g^{t2 + h)g<{h)g>{-h) 



X ^^(t2)sgn(t3) 



(A.15) 



dti 



dU 



X 



g<{-ti)g>{ti-t2-h)g>{U-t2) 
-9^{-ti)g''{ti-t2-h)g<{ti~h) 

g''{t2 + t3)g>{h)g<{-t^) 
-g^{h + h)g<{h)g>{-h) 



(A. 16) 



Next, the terms for diagrams denoted in Figs. 3.3 (j) and (k) are equivalent 
except for the spin indices and written by 



dti 



dU 



dU e 



iEti 



g>{ti)g<{-U)g>{U-h-h) 
-g^{h)g>{-h)g<{h-t2-h) 

g^{h)g>{h)g<{-h) 

xg^it2) +g''ih)g^ih)g>{-t3) 
+gHh)g^{t3)gH-h) 



(A.17) 



.a(4) 



X 



dti 



dU 



dU e 



iEti 



g<{h)g>{-h)g<{h-t2-h) 
-g^{h)g<{-ti)g>{h-t2-h) 

g^ih)g>its)g<i-h) 
+gHt3)gHt3)g''{-h) 
+gHt3)g^{t3)gH-t3) 



(A.I8) 



In addition, the terms for diagram illustrated in Fig. 3.3(1) are expressed by 



.r(4) 



X 



dti 



dto 



dt-i e 



iEti 



g>{ti)g>{ti)g<{-ti+t2 + t^) 
-g<(ti)g<{ti)g>{-ti+t2 + t^) 

^^<^h)g>{-t,)g<{t,) 



^g 



-t2) 



g 

^g' 
i-g' 



-t3)g^{-t3)g^{h) 
-t3)g^{-t3)g^{h) 



(A.19) 



.a(4) 



dti 



oo roo 

dh / 

-oo J —oo 



dU e*^*i 



X g 



g<{h)g<{h)g>{-h+h + h) 
-g>{ti)g>{ti)g<{-ti+h + h) 

gH-t3)g^i-t3)9^ih) 

+9^i-t3)9^{-t3)9^{t3) 



(A.20) 



6.2 Appendix B 

Expressions for Magnetization and Sus- 
ceptibility 

In the connection with the nonequihbrium state, the expressions for mag- 
netization and susceptibihty are written. Here, we consider the system where 
the magnetic field is apphed to the impurity. There the Zeeman term of the 
impurity, —BSz ( B is magnetic field ) is added to the Anderson Hamilto- 
nian. Magnetization for spin 1/2 is written by 

M = (Sz) = - {n,,)) = ^J dE[G<{E) - G<{E)l (B.l) 

where 

{h,^) = ^J dEG^E). 

For simplification, it is assumed that the system is noninteracting ( f/ = 
) and has symmetries: = Tr, {T^ + Tfi)/2 = T, and the applied voltage: 
I^L = eV/2, = -V/2. 

The Fermi distribution function can be written by 

/(x) = (e^ + l)-i = ^ + ^tangz). (B.2) 



From the formula of digamma function ip, 



1 

27ri 



1 X 
- + t — 

2 27r. 



1 . X 

2 ~^2^ 



- tan —t 



(B.3) 



Accordingly, the Fermi distribution function can be written in terms of 
digamma function by 



fix) = -- — {^1; 



1 X 
- + i — 

2 27r, 



1 .a; 

2 ~ ^2^, 



(B.4) 



The charge is expressed in terms of the Fermi distribution function using the 
residue theorem for finite magnetic field by 



{ndu) = f 



±B-ir 



(B.5) 



where T is temperature. If the right-hand side of Eq. ( B.5 ) is replaced with 
Eq. ( B.4 ), then, 



1 1 
1 1 



TT 



Imip 



1 ±B-iT 
- + i 

2 2nT 

±B + iV 



±B -iV 



2tiT 



2TiiT 



Magnetization at equilibrium state ( eV = ), therefore, reduces to 



M{B) = Imip 



TT 



1 B + iT 

2 ^ 2TiiT 



H — Im-^ 

TT 



-B + iT 



2mT 



(B.6) 



(B.7) 



For nonequilibrium state ( eVy^O ), the effective Fermi distribution function 
is obtained from Eq. ( 3.8 ) by 



ME) 



hiE) + fniE) 



using this, then, magnetization at eV^O is written by 
1 



M(5,e\/) = -| - -Im^; 

1 

-Imip 



+ 



1 B + eV/2 + tV 

2 + 

1 

2 ^ 27iiT 



2TiiT 
B + eV/2 + iV 



— iTSlIp 
TT 



H — Imi/j 

TT 



(B.8) 



1 B - eV/2 + iT' 

2 ^ 2TTiT 

1 -B - eV/2 + iT 

2 ^ 27riT 



The expression for nonequilibrium state is written as the sum of term of each 
chemical potential in the leads. 



When the limit is taken, the expressions are simplified as follows: in zero 
temperature limit, the expressions for magnetization M and susceptibility x 
at equilibrium state reduce to 



M{B) = - 

ZTT 



, B 

arctan ( — 



arctan 



-B 



X{B) 



1 [B 
—arctan — 

TT vr 

dM{B) 1 



dB 7r52 + r2. 
For nonequilibrium state( eV^Q ), 



M{B,eV) = — 
/vr 



(B + eV/2\ (B-eV/2' 
arctan + arctan 



(B.IO) 
(B.ll) 



[B.12) 



x{B,eV) 



V [S2 + {eV/2f + 



IX [{B + eV/2y + r2] [{B - eV/2f + 



[B.13) 



In isolated limit T— >-0, that the connection of the quantum dot with leads 
vanishes, the expressions for magnetization for thermal equilibrium reduce 
to 



M(5,T) 



1 


1 


2 


[e-^l^ + 


1 


/ B 


-tanh — 


2 


\2T 



[B.14) 



the Brillouin function as is known. In addition, susceptibility is also obtained 
by 



X(5,r) = — sech^ 



2T) 



[B.15) 



At nonequilibrium state, 



M{B,T,eV) = ^ 



tanh 



'B + eV/2' 
2T 



tanh 



'B-eV/2' 
2T 



(B.16) 



x{B,T,eV) 



8T 



2T 



2T 



In consequence, the expressions at nonequihbrium state are gained as the 
sum of term of each chemical potentiaL 
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